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Abstract: The classic Coupons-Collector Problem (CCP) is generalized. Only basic 
probability theory is used. Centerpiece rather is an algorithm that efficiently counts 
all fe-element transversals of a set system. 
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(N 
^ ■ 1 Introduction 

C^^ ■ Our setting is as follows. Let Vl^ be a set whose w many elements c will be viewed as "coupons" 

Let Ti = {Hi,- ■ ■ ,Hh} be any family of nonempty (not necessarily distinct) subsets. Thus 
[JJi Q W. By definition Hi contains exactly the coupons c of the i-th purpose. Put another 
way, each fixed coupon c € VF is multipurpose in the sense that it can serve many purposes i,j, 
and so forth, according to the sets Hi,Hj, ■ ■ ■ that contain c. If[j'H = W, then every coupon 
has at least one purpose. It is convenient to imagine the w many coupons as being located in 
an urn. 






In a length n trial a set of n coupons is picked at random one by one and all occuring purposes 
CN ■ are recorded. Thus for any picked coupon some of its purposes may have occured already and 

are not again taken into account. In a trial with replacement each coupon is put back into the 
f^ \ urn after its purposes have been ticked off. In a trial without replacement no drawn coupons are 

'nT ' put back (then necessarily n < w). A trial is successful if all h purposes show up. It is handy to 

call a successful trial sharply successful if the h purposes are only completed in the last draw. 

The Generalized Coupon- Collector Problem (GCCP) is to calculate the expected length £ of a 
sharply successful trial. We shall use the notation £ = lr{}i) and I = InriT'i) for a GCCP with 
replacement, respectively with no replacement. A prototypical example that can be cast as a 
GCCP with replacement, is the game of Roulette in section 6. 



% 



What is the classic Coupons- Collector Problem (CCP) and how does it fit our framework? The 
CCP can be phrased as follows: A firm encloses in each package of a certain product one coupon 
which is chosen at random from a stack of h different items. A customer who collects these 
coupons may gather some of them several times before he completes the whole set. Given the 
probabilities pi of picking a coupon of the ith type (or purpose) in a package, what then is the 
expected number (.{pi, ■ ■ ■ ,ph) of packages to be bought before all h types have occured? It is 
known that 

In particular, ii pi = P2 = ■ ■ ■ = Ph = j^ then (1) can be shown to simplify to 

^(i, •••,!) = /iHarm(/i) (2) 



where Harm(/i) := 1 + 2 + ' ' ' + Ti ^^ ^^^ harmonic number. 

For instance, setting /i = 6 in (2) one finds that a dice has to be thrown 14.7 times on average 
until all numbers have shown up. Applications of the CCP are plentiful, some of them outlined 
(with references) in [BH, section 1]. 

If without relevant loss of generality we assume that allpj's are rational numbers, we can build (in 

I J-f-\ 
various ways) a partitioning Ti* = {Hi, • • • , Hh} of a large enough u'-set W such that pi = — — 

w 
for all 1 < i < h. It follows that the CCP is a special case of GCCP and that 

e{pi,---,Ph) = ir{n*). (3) 

The £nr-analogue of CCP is of lesser interest. Observe that whereas (3) is well defined indepen- 
dent of Ti* (as long as pi = -^), this is not the case for lnr{T~i*)- For instance, returning to the 
dice problem: 

£.({{!}, {2}, •••,{6}}) = 4({{1,2},{3,4},...,{11,12}}) = 14.7 

1979 



.({{1},{2},...,{6}}) = 6, 4r({{l,2},{3,4},...,{ll,12}}) 



231 



Surprisingly perhaps, our approach to GCCP appeals more to the GCCP with no replacement 
(section 3). Only afterwards in section 4 we also tackle the GCCP with replacement. A numerical 
evaluation of our methods (in particular when applied to (1)) follows in section 6. Section 7 
tackles another problem treated in [BH], but again generalizes the perspective from coupons to 
purposes. Sections 2 and 5 handle the core of it all, i.e. is an efficient algorithm to count all 
fixed-cardinality transversals of a finite set system. 



2 It's all about counting transversals 



In this and the next section all trials are silently assumed to be with no replacement. Mathemat- 
ically our approach to the GCCP is straightforward; the challenge is its algorithmic realization in 
the next section. As to the mathematics, for fixed k G {0, 1, • • • , w} let qt be the probability that 
a length k trial is successful (in particular go = and q^ = 1). Observe that a successful trial 
induces a transversal X oil-Lva the sense that X r\Hj ^ $ for all 1 < j < h. In fact, if r^ is the 
number of /c-element transversals of H then exactly kWk trials among the w{w — l)---{w — k + l) 
many length k trials are successful, and so 

w[w — I) ■ ■ ■ [w — k + 1) 

Hence 

Sk ■■= Qk - qk~i {I <k <w) 

is the probability that a length k trial is sharply successful. Therefore 

inrCH) = lsi + 2s2H \- ws^, {= wq^ - qw-1 q2-qi) (5) 




3 The GCCP without replacement - inri'H) as a fraction 



By way of example we sketch an efficient algorithm for calculating (4) and hence (5). Consider a 
set W = {ci, • • • , cg} of eight coupons, each one of which serves between one and three purposes, 
according to Table 1: 
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Table 1 

For instance the trials ci, 03,05 and C6,C2, 03,07 are successful. The first is sharply successful, 
the second not. In order to calculate the expected length of a sharply successful trial, we read 
Table 1 row-wise and put H := {Hi,H2,H^,H^} with 



-H"! = {Cl, 02,03}, F2 = {03,06,07,08}, -^3 = {C2,C4, 05,06}, -^4 = {ci, O3, O4, 06, Os} 



As seen in Section 2 we need to count the r^ many fc-element transversals of Ti (1 < k < 
This is achieved by the transversal e-algorithm [W] which delivers the set of "models" 

Mod := {X QW : X is transversal of Ti} 

as a disjoint union of {0, 1, 2, e}- valued rows rf. 
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An idea of how these kind of {0, 1, 2, e}-valued rows arise will be given in section 5.1. Each row 
encodes a set of 0, 1-bitstrings of length 8 that correspond to subsets X QW contained in Mod. 
Here 2 is a "don't care" symbol, i.e. the corresponding entry can be or 1. Thus row r2 encodes 
2^ bitstrings; one of them is 10010100 which matches the transversals {01,04,06}. A string of 
symbols ee • • • e by definition means that any 0, 1-pattern with at least one 1 is allowed. In 
other words, only 00 • • • is forbidden. If serveral such e-patterns occur within one row, they are 



mutually independent and notationally distinguished. Thus row r^ contains (2 
bitstrings, one of them is 10011001. It follows that 



1)(2^ 



1) 
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|Mod| = 120 + 16 + 48 + 6 + 9 = 199 
For any {0, 1, 2, e}-valued row r on a set W and any nonnegative integer k G {1, 2, • • • , w} put 



Card(r, k) 



\{X er: \X\ = k}\. 



(6) 



If the set Mod of all transversals of some hypergraph 7^ on VF is represented as disjoint union 
of {0, 1, 2, e}-valued rows n, ^2, • ' ' > '<'R-, then 



Tfc = Card(ri,/c)+ Card(r2, fe) + • • • + Card(r/j, /c). 
In our toy example the numbers Card(rj, k) can be determined ad hoc and we get: 
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Tk = 
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Having the r^'s we can evaluate (4) and get: 
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Hence (4) gives 



^nrin) 



8qs - 97 



92 -91 



449 
140 



m = 3.207 



(8) 



In section 5 we will see how the above ad hoc manner gives way to a systematic procedure to 
calculate the numbers Card(r, A;) and thus (.nr{T~L)- 



4 The GCCP with replacement - ir{^) as a priori infinite sum 



Without further mention all trials in section 4 are with replacement. Analogous to section 2 we 
let t'^ be the number of successful length n trials, i.e. trials where all purposes of coupons have 
occured at some point (so t'^ = 0). Therefore, 



in 



4 (^ > 0) 



is the probability that a length n trial is successful, and 

S'n ■■= <i'n - 9^-1 {n > 1) 



(9) 



(10) 



is the probability that a length n trial is sharply successful. Hence the expected length of a 
sharply successful trial is 



<n) 



^ns'^ 



(11) 



n=l 



As to calculating the numbers t^, observe that no matter how coupons Cj are repeated in a 
length n trial, the underlying set of (distinct) coupons must be a /c-element transversal X of 7i, 
for some k < n. For a fixed set of coupons X = {ci, • • • , c^} C W the number of length n trials 
with underlying set X equals the number of ways to distribute n distinct balls (corresponding 
to the positions in the trial) to k distinct buckets (corresponding to the coupons) in such a way 
that no bucket stays empty. It is well known that this number is k\S{n, k) where S{n, k) is the 
Stirling number of the second kind. Accordingly, with Tj as in section 2, we deduce that 



t' 



' Ti{ll)S{n, 1) + T2{2l)S{n, 2) + • • • + Tn{n\)S{n, n), n<w 

(12) 
^ n{ll)S{n,l)+T2{2l)S{n,2) + --- + T^{wl)S{n,w), n>w 



q'u = ^ = 0.993319 



^1 
81 

Fortunately the infinite sum in (11) can be evaluated as a finite sum. What's more, one can 
express it in terms of the probabilities qk that a length /c-trial with no replacement is successful. 
Recall how q^ is coupled to r^ according to (4). 




For instance, for our example the Theorem yields ^r(^) = fl which is about 0.726 larger than 
^nri'H) from (8). Notice that d.ri'H) = wYiwva.{w) if and only if all qt = 0, which is the classic 
case where each coupon has only one purpose. The other extreme iri'H) = 1 occurs if and only 
if all gfc = 1, which means that every coupon fulfils every purpose. 

The Theorem is due to Svante Janson and (independently, with a more elaborate argument) to 
Stephen Wagner. Also Dirk Laurie contributed to the the matter. The next version of this draft 
will contain a full proof and hopefully feature some of the named people as co-authors. 



Calculating the number of /c-element transversals of a set sys- 
tem 



As mentioned article [W] features the fine details of the transversal e-algorithm, along with the 
calculation of the numbers Card(r, k) of any given {0, 1, 2, e}-valued row r. Rather than simply 
repeating this here, we take the opportunity to slightly extend transversals to transversouls. 
Correspondingly the transversal e-algorithm sketched in Section 3 becomes the transversoul e- 
algorithm. The gist of its row-splitting mechanism will be conveyed in section 5.1. The emphasis 
will however be (section 5.2) on how Card(r, k) can be computed for any single row r that the 
transversoul e-algorithm delivers. This works pretty much the same way as for transversals and 
the extra generality will enable us to tackle another coupons-related problem in section 7. 



5.1 The transversoul e-algorithm 

So let % = iHi, • • • , Hh} be a set system on 1^ = [w] and let a = (ai, • • • , ah) a fixed vector 
with integecl components aj > 1. If X <^W is such that 

{\fl<i<h) \XnHi\>ai (13) 

then X is called, tongue in cheek, a (positivqil) a-transversoul because it has more "soul" than 
an ordinary transversal where all ckj = 1. 

In order to present the transversoul e-algorithm in a nutshell, take W = [12] and let 

n := {{1,2, 3, 4, 5}, {6, 7, 8}, {4, 5, 6, 9, 10, 11, 12}} 
a := (ai,a2,a3) = (2,1,3) 

Akin to section 3 we wish to compactly encode the family Mod of all a-transversouls as a disjoint 
union of multivalued rows. The key idea is the symbolism e(s)e(s) • • • e(s), where s > 1 is an 
integer, and where the number m of symbols e(s) occuring must be greater than s. By definition 
only 0, 1-bitstrings X are allowed that have at least s entries 1 on the positions occupied by the 
symbols e(s). In particular (say) e(l)e(l)e(l) amounts the previously introduced eee. Because 

m 

^m — s + 1 
many old type e-constraints, for instance 



a constraint e(,s) ■ ■ ■ e(,s) of length m can always be written as an intersection of 



(e(2), e(2), e(2), e(2)) = (e, e, e, 2) n (e, e, 2, e) n (e, 2, e, e) n (2, e, e, e), 

employing e(s) • • • e(s) is mainly justified by reasons of succinctness and efficiency. As we shall 
see in a moment the use of "bubbles" e{s) ■ ■ ■ e{s) naturally brings about the use of bubbles 
g(s) ■ ■ ■ g{s) which by definition means that only bitstrings are allowed that have exactly s entries 
1 on the positions occupied by the symbols g{s). 

Coming back to our example, it is clear that if Ti was {Hi,H2} then the family of all (2,1)- 
transversouls X C W could be written as the row Fq in Table 1. (For systematic reasons we 
always write e(l) • • • e(l) rather than e • • • e.) In order to sieve exactly those X GTq that happen 
to be (2, 1, 3)-transversouls of "H we focus on HiCiH^ = {4, 5} and i?2nif3 = {6}. Combining the 
three options 00, ^(1)^(1), 11 of the former intersection with the two options 0, 1 of the latter 
intersection results in the six mutually disjoint rows ri to Tq. Consider e.g. row f4. Switching 
the last two components of e(2)e(2)e(2)e(2)e(2) in Yq to g{l)g{l) in ¥4^ forces its first three 
components to be e(l)e(l)e(l) in ¥4^. Similarly, switching the first component of e(l)e(l)e(l) in 
Tq to 1 in ¥4 frees its last two components to be 22 in ¥4. What's more, since g{l)g{l)l in ¥4 
means that either {4, 6} C X or {5, 6} C X for all X G ¥4, one needs to write e'(l)e'(l)e'(l)e'(l) 
at the end of ¥4 in order to ensure that \X D H^l > 3 for all X (^¥4. 



*One could allow ai — but that just amounts to dropping Hi from "H. 

^Here "positive" refers to the fact that all inequalities are of type >. This has the effect that the set Mod 
of all Q-transversouls is a set filter, and hence can be viewed as the set of all models of some positive Boolean 
function. Other ways to twist transversals to transversouls will be considered elsewhere. 
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Table 2 

As is the case for transversals [W, Thni.l] one can show that for given 7i and a also the N many 
a-transversouls X can be generated in time 0{Nw'^h?'). The same bound holds if merely the 
X's with \X\ > k need to be generated. The proof cannot be saved when \X\ = k is required 
(as it is in many applications) but in practise the transversoul e-algorithm performs as well as 
the transversal e-algorithm for which extensive experimental results are given in [W]. 



5.2 Focusing on individual rows 



It is an easy matter to calculate the (total) cardinality of a {0, 1, 2, e(,s)}-valued row; say 



Fsl 



(2^ - 1) . 2 • (2^ - 1) . (2^ - 5) 



462. 



Trickier is it to get Card(r, k) (defined as in (5)) for aW 1 < k < w. We proceed recursively. 
Assume that 



{e{b),---,e{b), 



,g{d),---g{d), e{s) , ■ ■ ■ , e{s)) 



is a typical row generated by the transversoul e-algorithm. Hence for some fixed a = (ai, ■ ■ ■ , ah) 
each X € r is an a-transversoul. Let 



(e(6),---,e(6), 



,g{d),--- ,g{d)) 



be the row obtained from r by dropping the last bubble, in our case an e-bubble. By induction 
all numbers 

card*(i) := |{X € r* : |X| = i}| (j < u;) 

are known. Thus card*(j) = if j is larger than the length of r*, and we put card*(j) := for 
j < 0. Letting m > s he the length of e(s) • • • e(s) one verifies that 



Card(r, k) 



^("']card*(fc-i) 

j=s ^ -^ ^ 



(14) 



If the last bubble of r was g{s) ■ ■ ■ g{s) rather than e(s) • • • e(s), then 

Card(r, fc) = (™) caid* (k — s) 



(15) 



As in [W, Thm.4] one shows that calculating Card(r, 1), Card(r, 2), • • • , Card(r, i^) this way 
costs 0{Kw^). 



6 Experimental results 



In this section we are only dealing with the transversal e-algorithm of Section 3, not its extension 
to transversouls. In subsection 6.1 it suffices to focus on just one {0, 1, 2, e}-valued row, but 
subsections 6.2 and 6.3 require a proper run of the transversal e-algorithm. 



6.1 Running the e-algorithm on the classic CCP 



Boneh and Hofri [BH, p. 43] emphasize the computational difficulty to evaluate (1) as h increases, 
and then go on to use integration for approximation. Recall that for rational pj's, say 

12 3 4 

P' = W'= 10' ^^ = W'= 10' 

our approach uses W = [10] and 

n = {{1},{2,3},{4,5,6},{7,8,9,10}}. 
Because the sets in Ti are disjoint, we can do with a single {0, 1, 2, e}-valued row 

r = (1, 62,62, 63,63,63 64,64,64,64) 

Using the technique of section 5.2 one readily computes Tk = Card(r, k) (1 < k < 14), and from 
them ir{T~i) according to the Theorem. 

Table 3 compares the 6-algorithm with the inclusion-exclusion approach (1) on such instances 
(pi, • • • ,ph) of the particular but natural type 
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which is uniquely defined by h (= first column in Table 3). As to inclusion-exclusion, we used 
a so called Gray-code in order to economically generate the subsets of [h] one by one from their 
predecessors, and also used that for common denominator probabilities one can simplify the 
terms in (1); say 
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Table 3 



The value of iri'H) is rounded to 6 digits albeit Mathematica, provided with the numbers 
Tk {^ ^ k < w), delivered the exact value as a fraction of two huge integers. For instance 



h = 150 gives w = 11325 and 24.5 sec of the 2508.2 sec total time were spent on plugging 
'Ti,T2,- ■ ■ T11325 into the Janson- Wagner formula. As is apparent, inclusion-exclusion (formula 
(1)) cannot compete. 



6.2 Roulette and some random instances 

In the popular game of Roulette a small metal bullet is spun and stopped at random on one 
of the w = 37 numbers 0, 1, 2, • • • , 36. Apart from each one of these numbers (coupons) has 
several properties (purposes) which head the columns of Table 4. As seen from the columns 
themselves, e.g. Manque and Passe (French) means the "first 18" and the "last 18" numbers. 

Thus the columns define the sets Hi {1 < i < h = 12) in H. Using the method of section 4 we 
found that the expected number of trials necessary to get all purposes is ^r(^) = 7.217 which 
is likely the correct rounded value (the cut off size in (11) was 0.00001). If after each trial one 
prevents the spinning wheel from delivering the same number again, then by (5) the expected 
number of trials is exactly inri'H) = 6.56233559797182. 
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Table 4 

Leaving Roulette, we also created (Table 5) some random instances with various parameters 
w, h, d. The larger the "density" d € (0, 1) the more wide spread the purposes are. For simplicity 
all Hi were chosen of equal cardinality \Hi\ = dw. 
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w h d 


4r(^) 


4(^) 


time 


10 5 0.3 


4.924 


6.6905 


0.015 


10 50 0.3 


7.127 


11.548 


0.094 


10 50 0.6 


3.992 


4.826 


0.14 


40 100 0.05 


36.534 


96.84 


29.9 


40 100 0.5 


7.109 


7.758 


716.9 


30 3500 0.1 


27.556 


70.485 


449.1 


20 3500 0.3 


13.4747 


21.524 


4619.9 



Table 5 



6.3 Information spreading and the expected time to conquer a chess board 

In many GCCP applications the purposes of a coupon c are other coupons, namely those that c 
wishes to "influence" in some way. More succinctly, we may consider a graph G whose vertex 
set W we imagine as a group of h = w people whose friendship relations are reflected by the 
edges of G. Suppose members c G W are phoned at random from outside W and told a piece 
of information. If c shares the news with all his friends, what is the expected number lriT~i) of 
phone calls necessary before the whole of W is informed? What is the analogue number inri'H) 
when nobody is phoned twice? 

A pleasant instance of the graph framework, where friendship turns to aggression though, is 
the problem to determine the average number of queens it takes when they are placed on a 
chessboard at random until all 64 squares are threatened. What are the corresponding numbers 
when only kings, or other figures (even of mixed types) are used? 



7 The likelihood of getting the z-th purpose at least ai times in 
k drawings 



For mere coupons instead of the more flexible purposes, and for trials with replacement, the 
question in the title is posed in [BH, p. 47] and a complex analysis involving the Cauchy integral 
formula ensues. Here we address the problem in a completely different manner. In particular, 
we consider trials without replacement, though trials with replacement could be handled as in 
section 4 by introducing Stirling numbers. 

Consider a set M^ = {ci,C2, • • • ,012} of twelve coupons, each one of which having one or two 
purposes as follows: 





Cl 


C2 


C3 


ca 


C5 


C6 


Cl 


C8 


eg 


ClO 


Cll 


C12 


purpose 1 


X 


X 


X 


X 


X 
















purpose 2 












X 


X 


X 










purpose 3 








X 


X 


X 






X 


X 


X 


X 



Table 6 
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Fix k G [12]. If the coupons are drawn with equal probabihty j^ ^^'^ without replacement, what 
is the likelihood Qk that after exactly k drawings the purposes 1,2,3 have shown up at least 
2, 1, 3 times respectively? Ad hoc one verifies that say 

Qi = Q2 = 0, Qs > 0, Qg< 1, Qio = Qii = Qu = 1. 

As to why Qs > 0, look at C4, C5, ce- In order to get the precise values of the QkS observe that 
Ti in section 5 matches Table 6. Analogous to Section 2 define 

Tk := number of k-element (2, 1, 3) — transversouls of H. 

Then as in (3) one gets (here w = 12): 

w[w — 1) ■ ■ ■ [w — k + I) 

Furthermore, it is clear that 

Tfc = Card(ri, k) + Card(r2, k) + ■ ■ ■ + Card(r6, k) 

where r^ is from Table 2. We e.g. focus on r^ and apply the procedure of 5.2 to obtain Card(r3, k) 
for all k € [12]. Hence put 

pi = e(l)e(l)e(l) 

P2 = e{l)e{l)e{l)g{l)g{l)0 

P3 = e(l)e(l)e(l)5(l)5(l)0e'(l)e'(l) 

rs = e(l)e(l)e(l)5(l)5(l)0e'(l)e'(l)e(2)e(2)e(2)e(2) 

For instance (15) was used to get the numbers of p2 from the ones of pi: 





1 


2 


3 


4 


5 


6 


7 


8 


9 


10 


11 


12 


Pi 


3 


3 


1 





























P2 





6 


6 


2 


























P3 








12 


18 


10 


2 




















n 














72 


156 


144 


70 


18 


2 









Table 7 

As another illustration, applying (14) to the numbers of ps we get 

Card(r3, 7) = (f\ ■ 10 + Q ■ 18 + (fj ■ 12 



144 



Similarly one gets the data for the other rows rj. The numbers T^'s result as the column sums, 
and the sought probabilities by using (16). 
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1 


2 


3 


4 


5 


6 


7 


8 


9 


10 


11 


12 


ri 

















24 


26 


9 


1 











f2 














18 


54 


61 


33 


9 


1 








rs 














72 


156 


144 


70 


18 


2 








Ti 











24 


108 


212 


238 


166 


72 


18 


2 





r-, 











8 


40 


86 


104 


77 


35 


9 


1 





re 








1 


9 


36 


84 


126 


126 


84 


36 


9 


1 


n 








1 


41 


274 


616 


699 


481 


219 


66 


12 


1 


Qk 








0.005 


0.083 


0.346 


0.667 


0.883 


0.972 


0.995 


1 


1 


1 



Table 8 
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COUPON COLLECTING WITH OR WITHOUT REPLACEMENT, AND 
WITH MULTIPURPOSE COUPONS 

Marcel Wild and Svante Janson 



Abstract: The classic Coupon-Collector Problem (CCP) is generalized. Only basic 
probability theory is used. Centerpiece rather is an algorithm that efficiently counts 
all /c-element transversals of a set system. 

o 

(N 
^ ■ 1 Introduction 

C^^ ■ Our setting is as follows. Let Vl^ be a set whose w many elements c will be viewed as "coupons" 

Let Ti = {Hi,- ■ ■ ,Hh} be any family of nonempty (not necessarily distinct) subsets. Thus 
[JJi Q W. By definition Hi contains exactly the coupons c of the i-th purpose. Put another 
way, each fixed coupon c € VF is multipurpose in the sense that it can serve many purposes i,j, 
and so forth, according to the sets Hi,Hj, ■ ■ ■ that contain c. lf[j'H = W, then every coupon 
has at least one purpose. It is convenient to imagine the w many coupons as being located in 
an urn. 






In a length n trial a set of n coupons is picked at random one by one and all occuring purposes 
CN ■ are recorded. Thus for any picked coupon some of its purposes may have occured already and 

are not again taken into account. In a trial with replacement each coupon is put back into the 
f^ \ urn after its purposes have been ticked off. In a trial without replacement no drawn coupons are 

'nT ' put back (then necessarily n < w). A trial is successful if all h purposes show up. It is handy to 

call a successful trial sharply successful if the h purposes are only completed in the last draw. 

The Generalized Coupon- Collector Problem (GCCP) is to calculate the expected length ^ of a 
sharply successful trial. We shall use the notation £ = lr{}i) and I = InriT'i) for a GCCP with 
replacement, respectively with no replacement. A prototypical example that can be cast as a 
GCCP with replacement, is the game of Roulette in section 6. 



% 



What is the classic Coupon- Collector Problem (CCP) and how does it fit our framework? The 
CCP can be phrased as follows: A firm encloses in each package of a certain product one coupon 
which is chosen at random from a stack of h different items. A customer who collects these 
coupons may gather some of them several times before he completes the whole set. Given the 
probabilities pi of picking a coupon of the ith type (or purpose) in a package, what then is the 
expected number (.{pi, ■ ■ ■ ,ph) of packages to be bought before all h types have occured? It is 
known that 

In particular, ii pi = P2 = ■ ■ ■ = Ph = j^ then (1) can be shown to simplify to 

^(i, •••,!) = /iHarm(/i) (2) 



where Harm(/i) := 1 + 2 + ' ' ' + TT ^^ ^^^ harmonic number. For instance, setting /i = 6 in (2) 
one finds that a dice has to be thrown 14.7 times on average until all numbers have shown up. 
Applications of the CCP are plentiful, some of them outlined (with references) in [BH, section 

!]• 

If without relevant loss of generality we assume that allpj's are rational numbers, we can build (in 

various ways) a partitioning 7i* = {Hi, • • • , Hh} of a large enough w-set W such that pi = 

for all 1 < i < /i. It follows that the CCP is a special case of GCCP and that 

e{pi,---,Ph) = ir{n*). (3) 

The ^„r-aiialogue of CCP is of lesser interest. Observe that whereas (3) is well defined indepen- 
dent of T-L* (as long as pi = -^), this is not the case for inriH*). For instance, returning to the 
dice problem: 

^,({{1},{2},...,{6}}) = 4({{1,2},{3,4},...,{11,12}}) = 14.7 

1979 
4,({{1},{2},...,{6}}) = 6, 4r({{l,2},{3,4},...,{ll,12}}) : 



231 



Surprisingly perhaps, our approach to GCCP appeals more to the GCCP with no replacement 
(section 3). Only afterwards in section 4 we also tackle the GCCP with replacement. A numerical 
evaluation of our methods (in particular when applied to (1)) follows in section 6. Section 7 
tackles another problem treated in [BH], but again generalizes the perspective from coupons 
to purposes. Sections 2 and 5 handle the core of it all, i.e. an efficient algorithm to count all 
fixed-cardinality transversals of a finite set system. 

This article is mainly about algorithms which are due to the first author; the second author's 
contribution is the probabilistic proof in Section 4. 



2 It's all about counting transversals 



In this and the next section all trials are silently assumed to be with no replacement. Mathemat- 
ically our approach to the GCCP is straightforward; the challenge is its algorithmic realization in 
the next section. As to the mathematics, for fixed k £ {0, 1, • • • , w} let qk be the probability that 
a length k trial is successful (in particular go = and g^ = 1). Observe that a successful trial 
induces a transversal X of 7^ in the sense that X CiHj ^ $ for all 1 < j < /i. In fact, if r^ is the 
number of fc-element transversals of H then exactly kWk trials among the w{w — l)---{w — k + l) 
many length k trials are successful, and so 

^'^ " w{w-l)---{w-k + l) I = ^^ 1 ^"^^ 

Hence 

Sk ■■= Qk - Qk-i (l <k <w) 

is the probability that a length k trial is sharply successful. Therefore 




-m 



Isi + 2s2 H h ws^ 



WQu 



lw-1 



92 - qi) 



(5) 



3 The GCCP without replacement - £nr{^) as a fraction 



By way of example we sketch an efficient algorithm for calculating (4) and hence (5). Consider a 
set W = {ci, • • • , cg} of eight coupons, each one of which serves between one and three purposes, 
according to Table 1: 





Cl 


C2 


C3 


ca 


C5 


C6 


c? 


cs 


purpose 1 


X 


X 


X 












purpose 2 






X 






X 


X 


X 


purpose 3 




X 




X 


X 


X 






purpose 4 


X 




X 


X 




X 




X 



Table 1 

For instance the trials ci, 03,05 and C6,C2,C8,C7 are successful. The first is sharply successful, 
the second not. In order to calculate the expected length of a sharply successful trial, we read 
Table 1 row-wise and put T-L := {/fi, if 2,-^^3,-^4} with 



-ffl = {C1,C2,C3}, F2 = {C3,C6,C7,C8}, -f^3 = {C2, C4, C5, Ce}, -f^4 = {ci, C3, C4, Ce, Cs} 



As seen in Section 2 we need to count the r^ many fc-element transversals of "H (1 < A; < 
This is achieved by the transversal e-algorithm [W] which delivers the set of "models" 

Mod := {X (IW : X \s transversal of %} 

as a disjoint union of {0, 1, 2, e}- valued (final) rows r^. 





Cl 


C2 


C3 


Ca 


C5 


C6 


C7 


C8 






ri 


2 


e 


1 


e 


e 


e 


2 


2 


ln| 


= 120 


r2 


1 








2 


2 


1 


2 


2 


1^2 1 


= 16 


n 


2 


1 





2 


2 


e 


2 


e 


ksl 


= 48 


ta 


e 


1 





e 


2 





1 





1^4 1 


= 6 


r5 


1 








e 


e 





e' 


e' 


k5l 


= 9 



An idea of how these kind of {0, 1, 2, e}-valued rows arise will be given in section 5.1. Each row 
encodes a set of 0, 1-bitstrings of length 8 that correspond to subsets X QW contained in Mod. 
Here 2 is a "don't care" symbol, i.e. the corresponding entry can be or 1. Thus row r2 encodes 
2^ bitstrings; one of them is 10010100 which matches the transversals {ci,ca,cq}. A string of 
symbols ee • • • e by definition means that any 0, 1-pattern with at least one 1 is allowed. In 
other words, only 00 • • • is forbidden. If serveral such e-patterns occur within one row, they are 
mutually independent and notationally distinguished. Thus row rs contains (2^ — 1)(2^ — 1) = 9 
bitstrings, one of them is 10011001. It follows that 



I Mod I 



120 + 16 + 48 + 6 + 9 



199 



For any {0, 1, 2, e}-valued row r on a set W and any nonnegative integer k G {1, 2, • • • , w} put 

Card(r,fe) := \{X £ r : \X\ = k}\. (6) 

If the set Mod of all transversals of some hypergraph ?^ on M^ is represented as disjoint union 
of {0, 1, 2, e}-valued rows n, r2, • • • , vr, then 

Tfc = Card(ri,/c)+ Card(r2, /c) + • • • + Card(r/j, /c). (7) 

In our toy example the numbers Card(rj, k) can be determined ad hoc and we get: 



k = 


1 


2 


3 


4 


5 


6 


7 


8 




\n\ 


Card(ri, A;) = 





4 


18 


34 


35 


21 


7 


1 




120 


Card(r2, k) = 





1 


4 


6 


4 


1 










16 


Card(r3, k) = 





2 


9 


16 


14 


6 


1 







48 


Card(r4, A;) = 








2 


3 


1 













6 


Card(r5, A;) = 








4 


4 


1 













9 
























Tk = 





7 


37 


63 


55 


28 


8 


1 




199 



Having the r^'s we can evaluate (4) and get: 



Qi 


Q2 


qs 


94 


95 


96 


97 


98 





1 

4 


37 
56 


9 

in 


5b 
56 


1 


1 


1 



Hence (4) gives 



inrin) 



898 - 97 



92 -91 



449 

140 



m « 3.207 



(8) 



In section 5 we will see how the above ad hoc manner gives way to a systematic procedure to 
calculate the numbers Card(r, k) and thus inriT~(-). 



4 The GCCP with replacement - ir{^) as a priori infinite sum 



Without further mention all trials in section 4 are with replacement. Analogous to section 2 we 
let t'^ be the number of successful length n trials, i.e. trials where all purposes of coupons have 
occured at some point (so t'^ = 0). Therefore, 



Jn 



-^ (n > 0) 

yjn \ - / 



is the probability that a length n trial is successful, and 

S'n ■■= q'n - Q'u-I (" > 1) 



(9) 



(10) 



is the probability that a length n trial is sharply successful. Hence the expected length of a 
sharply successful trial is 



\n) = Y.ns'^ (11) 



n=l 



As to calculating the numbers i^, observe that no matter how coupons Cj are repeated in a 
length n trial, the underlying set of (distinct) coupons must be a /c-element transversal X of 1-i, 
for some k < n. For a fixed set of coupons X = {ci, • • • , c^} C W the number of length n trials 
with underlying set X equals the number of ways to distribute n distinct balls (corresponding 
to the positions in the trial) to k distinct buckets (corresponding to the coupons) in such a way 
that no bucket stays empty. It is well known that this number is k\S{n, k) where S{n, k) is the 
Stirling number of the second kind. Accordingly, with Tj as in section 2, we deduce that 

r Tiil\)S{n,l)+T2i2\)S{n,2) + --- + Tnin\)Sin,n), n<w 

t'n = (12) 

[ Ti{ll)S{n,l)+T2{2l)S{n,2) + --- + TUwl)S{n,w), n>w 

Fortunately the infinite sum in (11) can be evaluated as a finite sum. What's more, one can 
express it in terms of the probabilities qk that a length A;-trial with no replacement is successful. 
Recall how q^ is coupled to r^ according to (4). 



(W~l \ UI — 1 

Harm(ti;) — \J = w}iaim.{w) — >^ /w-i\ • 

fe=i ^ / k=i V k ) 



Proof. Whether one arrives at a random sequence o" of A; distinct coupons by trials with or 
without replacement, in both cases a is sharply successful with the same probability s^ (as given 
in section 2). As opposed to the latter, in the relevant former case the expected number e^ of 

k 
steps required to get a is not k, but is well known [F, Example IX. 3(d)] to be e^ = > ■ 



w 



Recalling that (70 = and q^u = 1, we deduce: 

w w / k 



. ^ w + 1 — i 
1=1 



w 



w + 1 — i 

k=l k=l \ 1=1 / 

w fw w \ ( w w w 

[qi - qo)— + {q2 - gi) — H r H h{qw- qw~i) — H r H h — 

w \w w — lj \w w — 1 1 

—{q-w - qo) H z{qw -qi)^ \- T:\qw - qw-2) + -r{qw - q-w-i) 

w w — \ I 1 

UI— 1 



■u;Harm(it;) — wy 



qk 



w — k 

k=l 



By (4) we have r^ = qk{X} from which the rightmost formula in the Theorem follows. D 

For instance, for our example the Theorem yields i-r{T~i) = f| which is about 0.726 larger than 
^nri'H) from (8). Notice that ir{T~l') = wHaim{w) if and only if all qt = 0{k < w), which is 



the classic case where each coupon has only one purpose. The other extreme ^r(^) = w^ = 1 
occurs if and only if all qk = 1(1 < k < w) , which means that every coupon fulfils every purpose. 

Independently Stephan Wagner obtained the same result with a more complicated argument 

that tackled the Stirling numbers in (12) head on. Previously Dirk Laurie had initiated progress 

59 
in yet another way; he guessed the above — using Levin's f/-transformation and other tools 

15 
from numerical analysis. 



Calculating the number of /c-element transversals of a set sys- 
tem 



As mentioned article [W] features the fine details of the transversal e-algorithm, along with the 
calculation of the numbers Card(r, k) of any given {0, 1, 2, e}-valued row r. Rather than simply 
repeating this here, we take the opportunity to slightly extend transversals to transversouls. 
Correspondingly the transversal e-algorithm sketched in Section 3 becomes the transversoul e- 
algorithm. The gist of its row-splitting mechanism will be conveyed in section 5.1. The emphasis 
will however be (section 5.2) on how Card(r, k) can be computed for any single row r that the 
transversoul e-algorithm delivers. This works pretty much the same way as for transversals and 
the extra generality will enable us to tackle another coupons-related problem in section 7. 



5.1 The transversoul e-algorithm 

So let Ti = {Hi, ■ ■ ■ , Hh} be a set system oiiW = [w] and let a = (ai, • • • , ah) a fixed vector 
with integecl components a^ > 1. If X C VF is such that 

(Vl<i</i) \XnHi\>ai (13) 



then X is called, tongue in cheek, a (positivqil) a-transversoul because it has more "soul" than 
an ordinary transversal where all Oj = 1. 

In order to present the transversoul e-algorithm in a nutshell, take W = [12] and let 

n := {{1,2,3,4,5},{6,7,8}, {4,5,6,9,10,11,12}} 
a := (ai,a2,a3) = (2,1,3) 

Akin to section 3 we wish to compactly encode the family Mod of all o-transversouls as a disjoint 
union of multivalued rows. The key idea is the symbolism e(s)e(s) • • • e(s), where s > 1 is an 
integer, and where the number m of symbols e(s) occuring must be greater than s. By definition 
only 0, 1-bitstrings X are allowed that have at least s entries 1 on the positions occupied by the 
symbols e(s). In particular (say) e(l)e(l)e(l) amounts the previously introduced eee. Because 



*One could allow ai = but that just amounts to dropping Hi from "H. 

'Here "positive" refers to the fact that all inequalities are of type >. This has the effect that the set Mod 
of all Q-transversouls is a set filter, and hence can be viewed as the set of all models of some positive Boolean 
function. Other ways to twist transversals to transversouls will be considered elsewhere. 



m 



m 



s + 1 



a constraint e{s) ■ ■ ■ e{s) of length m can always be written as an intersection of 
many old type e-constraints, for instance 

(e(2), e(2), e(2), e(2)) = (e, e, e, 2) n (e, e, 2, e) n (e, 2, e, e) n (2, e, e, e), 

employing e(s) • • • e(s) is mainly justified by reasons of succinctness and efficiency. As we shall 
see in a moment the use of "bubbles" e{s) ■ ■ ■ e{s) naturally brings about the use of bubbles 
g{s) ■ ■ ■ g{s) which by definition means that only bitstrings are allowed that have exactly s entries 
1 on the positions occupied by the symbols g{s). 

Coming back to our example, it is clear that if H was {Hi,H2} then the family of all (2,1)- 
transversouls X Q W could be written as the row Fq in Table 1. (For systematic reasons we 
always write e(l) • • • e(l) rather than e • • • e.) In order to sieve exactly those X GYq that happen 
to be (2, 1, 3)-transversouls oiTi we focus on HiCiH^ = {4, 5} and H2r\H^ = {6}. Combining the 
three options 00, g{\)g{l)^ 11 of the former intersection with the two options 0, 1 of the latter 
intersection results in the six mutually disjoint rows ri to rg. Consider e.g. row r4. Switching 
the last two components of e(2)e(2)e(2)e(2)e(2) in Tq to g{l)g{\) in r4 forces its first three 
components to be e(l)e(l)e(l) in r4. Similarly, switching the first component of e(l)e(l)e(l) in 
Tq to 1 in r4 frees its last two components to be 22 in r^. What's more, since g{l)g{l)l in r4 
means that either {4, 6} C X or {5, 0} Q X for all X G r4, one needs to write e'(l)e'(l)e'(l)e'(l) 
at the end of T/^ in order to ensure that \X n H^\ > 3 for all X € r4. 
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2 


3 


4 


5 


6 


7 


8 


9 


10 


11 


12 


^'o = 


e(2) 


e(2) 


e(2) 


e(2) 


e(2) 


e(l) 


e(l) 


e(l) 


2 


2 


2 


2 




























ri = 


e(2) 


e(2) 


e(2) 











e(l) 


e(l) 


e(3) 


e(3) 


e(3) 


e(3) 


^2 = 


e(2) 


e(2) 


e(2) 








1 


2 


2 


e'(2) 


e'{2) 


e'(2) 


e'(2) 


^3 = 


e(l) 


e(l) 


e(l) 


g(l) 


g(l) 





e'(l) 


e'(l) 


e(2) 


e(2) 


e(2) 


e(2) 


ta = 


e(l) 


e(l) 


e(l) 


g(l) 


g(l) 


1 


2 


2 


e'(l) 


e'(l) 


e'(l) 


e'(l) 


n = 


2 


2 


2 


1 


1 





e(l) 


e(l) 


e'(l) 


e'(l) 


e'(l) 


e'(l) 


re = 


2 


2 


2 


1 


1 


1 


2 


2 


2 


2 


2 


2 



Table 2 

As is the case for transversals [W, Thm.l] one can show that for given 1-L and a also the N many 
a-transversouls X can be generated in time 0{N%iP'h'^). The same bound holds if merely the 
X's with \X\ > k need to be generated. The proof cannot be saved when \X\ = k is required 
(as it is in many applications) but in practise the transversoul e-algorithm performs as well as 
the transversal e-algorithm for which extensive experimental results are given in [W]. 



5.2 Focusing on individual rows 



It is an easy matter to calculate the (total) cardinality of a {0, 1, 2, e(s), g'(s)}-valued row; say 



F3 



(2^ - 1) • 2 • (2^ - 1) . (2^ - 5) 



462. 



The numbers Card(r, k) (defined as in (6)) are not much harder to obtain. It is easy to see that 
Card(r, k) is the coefficient at x^ of the polynomial pol(r, x) which is defined as the product 



of these factors: For each symbol 1 in r take a factor x, for each symbol 2 a factor 1 + x, for 
each constraint e(s) • • • e(s) of length m a factor (™)a;* + (s+i)^^'*^^ + • • • + (™)x'", and for each 
constraint g{s) ■ ■ ■ g{s) of length m a factor (J^^x'^. Thus, for instance 

pol(r3,x) = {3x + 3x'^ + x^){2x){2x + x'^){6x'^ + Ax^ + x"^) 
= 72x^ + 156x6 + 144x'^ + TOx^ + ISx^ + 2x^0 
and so e.g. Card(r3,7) = 144. 

6 Experimental results 



In this section we are only dealing with the transversal e-algorithm of Section 3, not its extension 
to transversouls. In subsection 6.1 it suffices to focus on just one {0, 1, 2, e}-valued row, but 
subsections 6.2 and 6.3 require a proper run of the transversal e-algorithm. 



6.1 Running the e-algorithm on the classic CCP 

Boneh and Hofri [BH, p. 43] emphasize the computational difficulty to evaluate (1) as h increases, 
and then go on to use integration for approximation. Recall that for rational Pi's, say 

12 3 4 

^' = To' ^^ = To' ^^ = To' ^^ = To' 

our approach uses W = [10] and 

n = {{1},{2,3},{4,5,6},{7,8,9,10}}. 
Because the sets in % are disjoint, we can do with a single {0, 1, 2, e}-valued row 

r = (1, 62,62, 63,63,63 64,64,64,64) 

Using the technique of section 5.2 one readily computes r^ = Card(r, k) {\ <k < 10), and from 
them iri'H) according to the Theorem. 

Table 3 compares the 6-algorithm with the inclusion-exclusion approach (1) on such instances 
(Pi; ■ ■ ■ ^Ph) of the particular but natural type 

1 2 /i / ^ , /i(/i + l) 

Pi = — , P2 = —, ■■■ ,Ph = — [ hence w = l^ h ft = 

WW w \ 2 

which is uniquely defined by h (= first column in Table 3). As to inclusion-exclusion, we used a 

so called Gray-code in order to more economically generate the subsets of [h] one by one from 

their predecessors, and also used that for common denominator probabilities one can simplify 

the terms in (1); say 

1 1 w 



Pi+Pj+Pk J: + ^ + ^ i + j + k' 



h 


tr{n) 


exclusion 


incl-excl. 


10 


68.9846 





0.2 


15 


150.606 





7.7 


27 


474.463 


0.3 


43193 


50 


1600.38 


4.1 


- 


100 


6338.75 


72 


- 


150 


14215.1 


455 


- 


200 


25229.5 


1829 


- 


400 


100667 


96272 


- 



Table 3 

The value of ^r('W) is rounded to 6 digits albeit Mathematica, provided with the numbers 
"Tk {^ ^ k < w), delivered the exact value as a fraction of two huge integers. For instance 
h = 400 gives w = 80200 and 3108 sec of the 96272 sec total time were spent on plugging 
Ti, T2, • • • T80200 i^to the formula of Section 4. As is apparent, inclusion-exclusion (formula (1)) 
cannot compete. 



6.2 Roulette and some random instances 



In the popular game of Roulette a small metal bullet is spun and stopped at random on one 
of the w = 2>7 numbers 0, 1, 2, • • • , 36. Apart from each one of these numbers (coupons) has 
several properties (purposes) which head the columns of Table 4. As seen from the columns 
themselves, e.g. Manque and Passe (French) means the "first 18" and the "last 18" numbers. 



Thus the columns define the sets Hi {1 < i < h = 12) in Ti. Using the Theorem of section 4 the 
expected number of trials necessary to get all purposes is exactly 

54728027202913 

; 7.201. 



7600186994400 



If after each trial one prevents the spinning wheel from delivering the same number again, then 
by (5) the expected number of trials is exactly 



L 



65774035502891 
10043104242600 



6.549. 





Pair 


Impair 


Rouge 


Noir 


Manque 


Passe 


p 


M 


D 


l.c 


2.C 


3.C 





























1 




X 


X 




X 




X 






X 






2 


X 






X 


X 




X 








X 




3 




X 


X 




X 




X 










X 


4 


X 






X 


X 




X 






X 






5 




X 


X 




X 




X 








X 




6 


X 






X 


X 




X 










X 


7 




X 


X 




X 




X 






X 






8 


X 






X 


X 




X 








X 




9 




X 


X 




X 




X 










X 


10 


X 






X 


X 




X 






X 






11 




X 




X 


X 




X 








X 




12 


X 




X 




X 




X 










X 


13 




X 




X 


X 






X 




X 






14 


X 




X 




X 






X 






X 




15 




X 




X 


X 






X 








X 


16 


X 




X 




X 






X 




X 






17 




X 




X 


X 






X 






X 




18 


X 




X 




X 






X 








X 


19 




X 




X 




X 




X 




X 






20 


X 






X 




X 




X 






X 




21 




X 


X 






X 




X 








X 


22 


X 






X 




X 




X 




X 






23 




X 


X 






X 




X 






X 




24 


X 






X 




X 




X 








X 


25 




X 


X 






X 






X 


X 






26 


X 






X 




X 






X 




X 




27 




X 


X 






X 






X 






X 


28 


X 




X 






X 






X 


X 






29 




X 




X 




X 






X 




X 




30 


X 




X 






X 






X 






X 


31 




X 




X 




X 






X 


X 






32 


X 




X 






X 






X 




X 




33 




X 




X 




X 






X 






X 


34 


X 




X 






X 






X 


X 






35 




X 




X 




X 






X 




X 




36 


X 




X 






X 






X 






X 



Table 4 

Leaving Roulette, we also created (Table 5) some random instances with various parameters 
w,h,d. For simplicity all Hi were chosen of equal cardinality \Hi\ = dw. Thus d = 0.3 means 
that each one of the w coupons serves a random 30% of the h purposes. As to "final rows", 
recall Section 3. 



10 



w 


h 


d 


enrin) 


eAn) 


time (sec) 


final rows 


40 


10 


0.3 


7.69 


8.57 


1.6 


579 


40 


10 


0.7 


2.78 


2.86 


0.3 


123 


40 


40 


0.3 


10.61 


12.28 


686.3 


274 117 


40 


40 


0.7 


3.85 


4.00 


6.0 


2092 


40 


100 


0.3 


12.34 


14.66 


18349.4 


5 429 569 


40 


100 


0.7 


4.46 


4.68 


44.5 


9842 


100 


100 


0.7 


4.66 


4.75 


2435.5 


361428 



Table 5 



6.3 Information spreading and the expected time to conquer a chess board 

In many GCCP applications the purposes of a coupon c are other coupons, namely those that c 
wishes to "influence" in some way. More succinctly, we may consider a graph G whose vertex 
set W we imagine as a group of h = w people whose friendship relations are reflected by the 
edges of G. Suppose members c G W are phoned at random from outside W and told a piece 
of information. If c shares the news with all his friends, what is the expected number iriT~(-) of 
phone calls necessary before the whole of W is informed? What is the analogue number i.nr{T~i) 
when nobody is phoned twice? 

A pleasant instance of the graph framework, where friendship turns to aggression though, is 
the problem to determine the average number of queens it takes when they are placed on a 
chessboard at random until all 64 squares are threatened. What are the corresponding numbers 
when only kings, or other figures (even of mixed types) are used? 



7 The likelihood of getting the z-th purpose at least ai times in 
k drawings 



For mere coupons instead of the more flexible purposes, and for trials with replacement, the 
question in the title is posed in [BH, p. 47] and a complex analysis involving the Cauchy integral 
formula ensues. Here we address the problem in a completely different manner. In particular, 
we consider trials without replacement, though trials with replacement could perhaps be handled 
similar to section 4. 

Consider a set W = {ci,C2, • • • ,012} of twelve coupons, each one of which having one or two 
purposes as follows: 





Cl 


C2 


C3 


ca 


C5 


C6 


C7 


eg 


eg 


ClO 


Cll 


C12 


purpose 1 


X 


X 


X 


X 


X 
















purpose 2 












X 


X 


X 










purpose 3 








X 


X 


X 






X 


X 


X 


X 



Table 6 
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Fix k G [12]. If the coupons are drawn with equal probabihty j^ ^^'^ without replacement, what 
is the likelihood Qk that after exactly k drawings the purposes 1,2,3 have shown up at least 
2, 1, 3 times respectively? Ad hoc one verifies that say 

Qi = Q2 = 0, Qs > 0, Qg< 1, Qio = Qii = Qu = 1. 

As to why Qs > 0, look at C4, C5, ce- In order to get the precise values of the QkS observe that 
Ti in section 5 matches Table 6. Analogous to Section 2 define 

Tk := number of k-element (2, 1, 3) — transversouls of H. 

Then as in (3) one gets (here w = 12): 

Q^ = ~f v\ — ? rVT{ (16) 

w[w — 1) ■ ■ ■ [w — k + I) 

Furthermore, it is clear that 

Tfc = Card(ri, k) + Card(r2, k) + ■ ■ ■ + Card(r6, k) 

where r^ is from Table 2. Thus if we first list in the ith. row (1 < i < 6) the coefficients of 
pol(ri,x) (see 5.2) then the numbers Tk result as column sums: 





1 


2 


3 


4 


5 


6 


7 


8 


9 


10 


11 


12 


ri 

















24 


26 


9 


1 











r2 














18 


54 


61 


33 


9 


1 








rs 














72 


156 


144 


70 


18 


2 








ta 











24 


108 


212 


238 


166 


72 


18 


2 





r5 











8 


40 


86 


104 


77 


35 


9 


1 





re 








1 


9 


36 


84 


126 


126 


84 


36 


9 


1 


Tk 








1 


41 


274 


616 


699 


481 


219 


66 


12 


1 


Qk 








0.005 


0.083 


0.346 


0.667 


0.883 


0.972 


0.995 


1 


1 


1 



Table 8 
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